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Our aim is studying the thermodynamics of cosmological models including initial and final de- 
Sitter eras. For this propose, bearing Cai-Kim temperature in mind, we investigate the thermo¬ 
dynamic properties of a dark energy candidate with variable energy density, and show that the 
state parameter of this dark energy candidate should obey the uod ^ —1 constraint, whiles there is 
no interaction between the fluids filled the universe, and the universe is not in the de-Sitter eras. 
Additionally, based on thermal fluctuation theory, we study the possibility of inducing fluctuations 
to the entropy of the dark energy candidate due to a mutual interaction between the cosmos sectors. 
Therefore, we find a relation between the thermal fluctuations and the mutual interaction between 
the cosmos sectors, whiles the dark energy candidate has a varying energy density. We point to 
models in which a gravitationally induced particle production process leads to change the expansion 
eras, whiles the corresponding pressure is considered as the cause of current acceleration phase. We 
study its thermodynamics, and show that such processes may also leave thermal fluctuations into 
the system. We also find an expression between the thermal fluctuations and the particle produc¬ 
tion rate. Finally, we use Hayward-Kodama temperature to get a relation for the horizon entropy 
in models including the gravitationally induced particle production process. Our study shows that 
the first law of thermodynamics is available on the apparent horizon whiles, the gravitationally 
induced particle production process, as the dark energy candidate, may add an additional term to 
the Bekenstein limit of the horizon. 


I. INTRODUCTION 


The universe expansion is modeled by the so-called FRW metric 


ds 2 = dt 2 


a 2 (t) 


dr 2 

1 — kr 2 



( 1 ) 


describing a homogenous and isotropic cosmos. In this metric, a(t) is the scale factor while k = —1, 0,1 is the curvature 
parameter corresponding to open, flat and closed universes respectively. WMAP data indicates a flat universe p]|. 
Conformal form of this metric can be used to model the inhomogeneity of cosmos due to the structure abundance in 
scales smaller than 100 Mpc The universe expansion, in the standard cosmology (AC DAI), consists of four eras: 

(i) the primary inflation era ( a(t ) ~ e Ht ) exactly began after the big bang needed to remove the horizon problem. 

(ii) the radiation dominated era ( a(t ) ~ tz) finished by leaving an attractive trace which is now called either CMB 
or LSS. (iii) the matter dominated era (a(t) ~ fs) which has the major part in the structure formation process, 
(iv) finally, in agreement with the current phase of the expansion, the universe is undergoing an accelerated phase 
(a(t) ~ e Ht ) which is similar to the primary inflation era, where H = ^ and dot are the Hubble parameter and 
derivative with respect to time, respectively. In order to explain both the primary inflation and the current phase 
of the expansion leading to the coincidence and cosmological constant problems [lj we need an abnormal matters. 
Thermodynamic considerations show that the universe maintains this current phase of the expansion forever leading 
to avoid the big crunch problem ||. We should note that this model (A CDM) is in a very good agreement with the 
observations [1], |_lj In . It should be noted that if one accepts the thermodynamic predictions about the current stage 
of the expansion 3, then three questions including the big bang, coincidence and cosmological constant problems 
remain unsolved. 

There are various models proposed to eliminate such weaknesses of the standard cosmology by introducing either a 
new degree of freedom or a new parameter leading to explanations for the generator of the current acceleration phase 
(called dark energy (DE)) in the Einstein relativity framework &G1- Moreover, it is shown that the DE candidates 
may add an additional term to the Bekenstein entropy of the apparent and trapping horizons of the FRW universe 
19-121]. Recently, two models are introduced in the literatures to solve the mentioned problems of the standard 
cosmology whenever their physics are completely different from those of previous works |22I.|23|]. In these models, the 
universe expansion is began from an unstable de-Sitter spacetime leading to solve the horizon problem. In addition 
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to this, the final stage of the universe expansion is an eternal de-Sitter phase which is due to the thermodynamic 
equilibrium conditions and leading to solve the big crunch singularity proble_m [23 . l25j . Because of different physics 
behind these models, they become interesting for further investigations Ml, As the main difference of these two 
models, while a varying vacuum plays the role of DE in the model proposed in 22], the gravitationally induced particle 
production process is the backbone of the second proposal [23|,[25|. Moreover, it is shown that the first model satisfies 
the Bekenstein limit of entropy [21] . 

It seems that the apparent horizon of the FRW universe, as a marginally trapped surface which is located at 


ta 


1 



( 2 ) 


can be considered as a causal boundary for this spacetime [35l - l40i ] . Surface gravity of the apparent horizon is evaluated 
as 


n = ^J— d a (\/^hh ab d b C), 


( 3 ) 


leading to 


H M H 
K ~ _ 27t 1 27J2 ’ 


( 4 ) 


and we get 


T = — = - — (l + -^o 
2 tt 2tt 2 H 2 


), 


( 5 ) 


as the temperature on the apparent horizon [331 . 134j . Indeed, this temperature is called Hayward-Kodama temperature 
used to show the consistency between the Friedmann equations and the first law of thermodynamics [iol - [43j . Moreover, 
some authors define T = — ■M-, called Cai-Kim temperature, in order to get the positive value for temperature 


44]. Another motivation for T = 21, signalling us that this temperature may be considered as the temperature for 


fields confined to the apparent horizon, can be found in ref. [45|. Finally, it is useful to note that these definitions of 
temperature could not attract a common agreement I4ll - l43l ]. 

Bearing the various definitions of temperature in mind, since the apparent horizon can be considered as the causal 
boundary, some authors have been shown that the validity of the first law of thermodynamics on the apparent horizon 
leads to Friedmann equations [40l-l5p . In addition, it seems that it is necessary to consider a DE |i0 or 
modifying the Einstein equations [551 ] in order to be compatible with recent observations imposing the a(t) > 0 and 
a(t) > 0 conditions on the scale factor [56l - l59[ . These data are also in agreement with the generalized second law of 
thermodynamics in numerous models of DE providing an eternal thermodynamic equilibrium state for the universe 
3.l2if60M62|. More studies on the thermodynamics of final state of the cosmos can bejound in ref. [63. f64j. 

Moreover, observations admit an interaction between the dark sectors of the cosmos [65l - l7ll ]. Additionally, it seems 
that the mutual interaction between the dark sectors of the universe may solve the coincidence problem [7l! - l78[ . 
Considering thermal fluctuation theory in mind [79] , authors have shown that the entropy of event horizon is modified 
by a logarithmic correction [80j . In the cosmological setup, by using the Cai-Kim temperature, it is shown that 
these fluctuations may be interpreted as the result of a mutual interaction between the dark components [78j . The 
generalization of this approach to more cosmological models can be found in [sT| - [ 8G| . Therefore, it seems that one 
can find an expression for the mutual interaction between the dark sectors of the cosmos by taking into account 
the thermal fluctuations of the universe components. In addition, it is shown that a mutual interaction between the 
DE candidate and the other parts of cosmos may add an additional term to the Bekenstein entropy of the apparent 
horizon [2lt |. 

Here, according to the foregoing discussion and by following the approach considered in refs. [78l. IsTl - fSbl . we focus 
on the models proposed in Refs. [2ll. l23j. and try to find the suitable thermodynamic interpretations for these models 
by using the Cai-Kim temperature as well as the thermal fluctuation theory. In fact, we try to get a thermodynamic 
interpretation for the mutual interaction between the DE candidate and the other parts of cosmos. We also use 
the Hayward-Kodama temperature, and show that the gravitationally induced particle production process, as the 
DE candidate, may add an additional term to the Bekenstein entropy of the apparent horizon. The latter is in 
agreement with the prev ious attempts which predict that the DE candidates may correct the Bekenstein bound in the 
cosmological setups Il9l - I21| . For simplicity, we take G = h = c = 1 throughout this paper and we restrict ourselves to 
the k = 0 case. Moreover, dot denotes derivative with respect to time. 
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The paper is organized as follows. In next section, we decompose the universe sectors into a varying DE candidate 
and other parts with total density p and state parameter u> whiles, the cosmos sectors do not interact with each 
other. Then, we study the thermodynamics of the DE candidate in this model. In section III, we consider the model 
introduced in (22j |. and give a thermodynamical interpretation for this model using thermal fluctuations theory. To 
do this, the model will be considered as a model in which the cosmos sectors are interacting with each other. In 
section IV, we focus on the model proposed in Ref. @1, and show that the gravitationally induced particle production 
process leads to the thermal fluctuations and finally, we get a thermodynamical description for this model. We also 
show that the gravitationally induced particle production process may add an additional term to the Bekenstein limit 
of the apparent horizon of the flat FRW universe in section V. Last section is devoted to summary and concluding 
remarks. 


II. THERMODYNAMICAL DESCRIPTION OF DE MODELS WITH NON-CONSTANT ENERGY 

DENSITY 


For the flat FRW universe supported a DE candidate, Friedmann equations lead to 

H 2 = P + Pd), 


( 6 ) 


and 


— 2- H 2 = 8nG(p + p D ) 

a 


( 7 ) 


where po and pn are the density of DE and its corresponding pressure, respectively. In addition, p comprises 
other parts of cosmos which may include radiation, pressureless matter, dark matter and etc. p is also the pressure 
corresponding to the density p. Consider a DE candidate with profile density 


Pd{H) = =a+P H 2 + 1 H 2n , 

87r 


( 8 ) 


which converges to that of the ghost dark energy model by substituting a = (3 = 0 and n = \ @. Moreover, it 
covers the profile density of generalized ghost dark energy model by inserting a = 0 together with n = 1 fl7l ]. The 
cosmological constant model, as the trivial limit, is obtainable by inserting /3 = 7 = 0 flj. More models in which 
authors used a dynamic DE model to explain the current expanding phase can be achieved by choosing proper values 
for n, a, j3 and 7 [22I (87l - [93| . It is also shown that this profile density may add an additional term to the Bekenstein 
entropy of the trapping and apparent horizons of the FRW universe fljMbll ] . Using this equation, one gets 


P=H 2 (—-/3)-a- 7 H 2n . 

8n 

For the total energy momentum tensor, the energy momentum conservation law implies 

p + pd + 3 H (p(l + w) + pd{ 1 + 3w_d)) = 0, 

where dot denotes derivative with respect to time. This equation can be decomposed into 


and 


p + 3-ff p(\ + w) — 0, 


Pd + 3IL/7d( 1 + wo) = 0, 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


where 07 = y and pi are the state parameter and the pressure of the i ^ 1 sector, respectively. Such well-advised 
decompositions are valid whenever there is no mutual interaction between the dark sectors. In such a situation, from 
Eqs. (El) and G2D, clearly we have pc ^ 0 yielding u>d 7 ^ —1- For models introduced in Ref. [22J], the energy profile 
density of varying vacuum, as the DE candidate, is also given by ©, whiles u>d = — 1 , and satisfies the thermodynamic 
equilibrium conditions in current acceleration phase of the universe expansion [2_4]. It is shown that such models may 
avoid the big bang and big crunch singularities as well as the horizon problem, and can also provide a complete 
description for the history expansion of the universe [22j ■ We should note here that since the state parameter of the 
DE candidate satisfies the ojd = — 1 constraint in model proposed by Lima et al. [22], the decomposition of uni) 
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into m and © is possible if p = 0 leading to H = 0 because p = jjjH. Briefly, since model proposed by Lima 
et al. [22j satisfies the thermodynamic equilibrium conditions in the current de-Sitter accelerating phase (H = 0) 
H|, the decomposition of ijTUjl into Eqs. m and (1121) as the result of marginally thermodynamic equilibrium states 
(de-Sitter spacetimes with H = 0) of the model is reasonable. 

Derivation from Eq. ([6]) with respect to t and using m, leads to (Raychaudhuri equation) 

H =^-=-4n[p(l + ui) +p D (l+u D )\- (13) 

Now, if we define p c = 3-^- and use Eqs. (0, (0 and (flUl) we obtain 


H = -4:Tt[(ujd - u>)(a + f3H 2 + yif 2 ") + ^-(1 + w))]. (14) 

o7T 

Since we have only used Eqs. m-m in order to derive Eqs. m and ED, we should note that these equations are 
independent of the validity of Eqs. ED and CEb and thus the probable mutual interaction between the dark sectors. 
Indeed, the validity of equation ED. and therefore ED, is due to the Bianchi identity or the conservation of the total 
energy momentum tensor ED- Finally, by inserting p c into the Friedmann equation we get 

1 = Q d + Q, (15) 

where Dj = is the fractional energy density of the i^ 1 component of the cosmos. For the DE candidate confined to 

Pc _ 

the flat FRW universe enclosed by the apparent horizon 0, the Gibb’s law implies 


TdS D = dE D +p D dV. (16) 

In this equation, Sd is associated entropy to the DE while V = ^-r\ and Ed = p^V are the volume of the flat 
FRW universe and the energy of DE, respectively. Additionally, thermodynamic equilibrium condition implies that 
T (the temperature of DE) should has the same value as the temperature of the apparent horizon. Moreover, since it 
is unreasonable to have a fluid with negative temperature, meaning that T > 0, the Cai-Kim temperature may be a 
good option for the temperature of the DE candidate [U [3§, SI, S3, SE ZSi I8ll-l86l] 



1 

2irrA ’ 


where ryi is the apparent horizon radius of the flat FRW universe 0 



Therefore, for the volume and energy differentials we reach 


(17) 


(18) 


dV = 4:Tr(fA) 2 drA = —47 tH 4 dH 1 


(19) 


and 


dE D = pndV + Vdpo , 


( 20 ) 


leading to 


dS D = 2Trf A {pD(l + UJ D nTr{f A ) 2 dfA + —{rA) 3 dp D ), 


where we have used pd = Pd^d- In addition, since 


dpD dr a dH 
^ D dr a dH dt 

we find 

dpD _ 3 / 9 d (1 + wd ) 

dr A 47rr^(p(l + w) + p D (l +w.d))’ 


( 21 ) 


( 22 ) 


(23) 
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where we have used Eqs. ini), da and (HI. By combining this equation with (1211) and after some algebra, we obtain 

1 


dS" D = 8n'p" D (l + ^ D my(l- 


4^) 2 (p°(l+cu°)+p° D (l+^)) 




(24) 


We must note that the superscript (0) is used to indicate that this result is valid whenever there is no interaction 
between the cosmos sectors of the universe and therefore, f° A = -W. Since we did not use Eq. © to obtain this 
equation, this result (Eq. (l24l) ) is also valid in every cosmological model with the same Friedmann equation as Eq. © . 
Using Eqs. ©, m and CHI the above equation can be rewritten as follows 


dS° D = - 


87r 2 (a + /3 -Hq 7 -ffo”)(l + u>p) 

rj5 


( 1 - 


H o 


4t t[{u° d - w°)(a + PHi + ytf 2 ") + ^(1 + < 


-)dH 0 . 


(25) 


Now, using Eqs. © and (fl5l) to get 


dS° D _ 3^(1+^) _2_ 

dH 0 Hi 4[(w° -w 0 )f2°j + (l+a; 0 )] h 


(26) 


where the superscript/subscript (0) implies non-interacting case. Therefore, we find an expression for the entropy of 
the DE candidate when its energy density profile is varying with time as introduced in Eq. ©. It is useful to be 
noted that, in the u> 0 D = — 1 limit, dSk = 0 is obtainable meaning that pu = 0 which can be considered as either the 

or the marginally thermodynamic equilibrium states of model proposed by 
0 _ -j dS° n d 2 S° r 


fl, 


whiles Up = — 1, 


cosmological constant model of DE 

Lima et al. [2J. In fact, from (l26l) it is apparent that, wim^ — d[[ 

the thermodynamic equilibrium conditions are marginally satisfied in these eras 9 

initial and final de-Sitter spacetimes of models proposed by Lima et al. p2| are marginally thermodynamic equilibrium 
states. In the next section, we show that whenever a mutual interaction between the cosmos sectors moves the cosmos 
between these two marginally thermodynamic equilibrium states 


= 0 is available, meaning that 
The latter signals us that the 


22 ], it leaves a thermal fluctuations into the system. 


III. THERMODYNAMIC DESCRIPTION OF THE INTERACTING DE MODELS WITH 

NON-CONSTANT ENERGY DENSITY 

Here, we study the thermodynamics of a universe filled by a varying vacuum, as the DE candidate, together with 
another source of energy with density p interacting with each other whenever the energy density of the DE candidate 
is the same as the previous section ©. In this manner, decomposition of the energy-momentum conservation law (1101) 
implies 


p + 3Hp(l+u) = -Q, (27) 

and 

Pd + 3Hpr>(l + u;d) = Q, (28) 

where Q is the mutual interaction between the cosmos sectors. Clearly, ujd = —1 is accessible in this model while 
Pd ik 0. Some authors have considered this possibility (cod = —1), and showed that the model leads to the compatible 
outcomes with observational data and the thermodynamic equilibrium conditions [22l . [24| . In addition, it is also shown 
that, independent of the profile density ofthe DE candidate, the apparent horizon satisfies the Bckenstein entropy, 
in the interacting models with uid = — 1 [2lj ]. It seems that the mutual interaction between the dark sectors induce 
some fluctuations into the thermodynamic properties of systems which can be investigated by thermal fluctuations 
theory {78i. [80| - l8fj . Due to these fluctuations, the entropy is changed to Sd and it can be expanded as [7814861 ] 

Sd = Sd + Sd + Sd (29) 

In this equation, S^ is the entropy of the DE candidate when there is no mutual interaction between the cosmos sectors 
of the universe and S'/, = — | In CTq is logarithmic correction to the entropy where C = To/jyP- is the dimensionless 
heat capacity. S/, also concerns higher order terms. It is also useful to mention that this analysis is valid for all 
thermodynamical systems [80j |. Combining Eqs. (Hi and we get 

r - r dSp pMI _ ?go n 

° dT 0 p°(l + w°) + p° D (l + uj° d ) 


(30) 
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From Eq. 03 we have also 


C = Tr 


dS° r , 


D 


dT 0 


= H ( 


dS° r 


D 


dH 0 


leading to 


^ 87r 2 (a + (3Hq + 7FTg ra )(l + u>° D ) n 

“ 773 

n 0 


H o 


3 H? 


- w°)(a + !3Hl + 7 Hg n ) + -^(l + w 0 )] 


(31) 


(32) 


where we have used Eq. (l25l) in order to evaluate the heat capacity of the DE candidate with profile density (0. 
Therefore, one gets 


d !, n ^\a + fiH 2 + 1 H 2n ){l+^ D ) / 
ln: 777 V 


H 0 


Hi 


Att[{u° d - w°)(a + /3H 2 + yi/ 2 ™) + ^(1 + w 0 )] 

In deriving this equation we used again Eqs. (El and El along as (l32l) . It is a matter of calculation to show that 
?! 


- 1)]). (33) 


dSjj _ 12pH 0 + 2n'yH 0 

dH 0 ~ 2 a + /3H 2 + 7 i7 2n 


2 "- 1 2ff 0 - 4 tt[(lo 0 d - ui°)(2/3H 0 + 2n 7 g 0 2ri - 1 ) + ^(1 + w °)] 


H o - 4tt[(u;° - w°)(a + /3ff 0 2 + 7^") + ^(1 + <*> 0 )] 


_ [(u° D ~ cc°)(2/3g 0 + 2 n^- 1 ) + ^(1 +a; 0 )] 2 . 

[K - w°)(a + /3H 2 + 7^") + ^(1 + <*> 0 )] ' 

which can be simplified, by using Eq. m , as 

3Hqi 


(34) 


dSb _ 1 A&y 2 Hq - faH - u>°)(p° D y + ffi(l + a; 0 )] _ [(o,g, - u; 0 ) W + 1+cc 0 )] _ _2_ 


dHn 


Pd Hi - 4tt[(u;° - o^pO, + ^(1 + u;°)] 




ffn 


(35) 


where the prime stands for the derivative with respect to H$. Again, we should note that the subscript/superscript (0) 
is used to emphasize the non-interacting parameters. By using Friedmann equation 0 along as Eqs. El and 
it is easy to show that 


dpD 


Q — 3Hpz>(l + wp) 


dr a 47rr^(p(l + w) + pzj(l + u> D )) ‘ 

By following the recipe of previous section we get 

Q — 3Hpr>(l + u>d) 


dS D = 87r 2 p£)(l + u] D )(f A ) (1 + 


12nr A p D (l + wd)(p(1 + w) + p D {\ + ui D )) 


)df A . 


Using Eq. (ITS1) . we find 

dS D _ 87t 2 pd(1 + u>d) 


dH 


H 5 


3H 2 p D (l+u D )-QH 


12np D (l + uj D )(p( 1 + u>) + p D ( 1 + ud)) 


-!)• 


(36) 


(37) 


(38) 


As a desired result, the results of previous section are obtainable by substituting Q = 0. Bearing the model investigated 
in [22], [24[ in mind, where = — 1, and use Eq. 0 to obtain 


dS D 

dH 


27 tQ 


3(1 + u)H 4 {a + 7 H 2n - H 2 (± - 0)) 


(39) 


Comparing this equation together with (1261) , we can conclude that the mutual interaction may change the entropy of the 
varying vacuum leading to separation from the marginally thermodynamic equilibrium situation ( dS^ = d 2 S^ = 0). 
Since Q = Pd in this model, by using 0 and El, we get 


3FP 

87T 


)• 


Q = 47 Tp' D (u + 1 )(pd 


(40) 
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Considering the definition of the fractional energy density, this equation can written as 

Q = + 1 )H 2 (VLd — 1 ), 


(41) 


leading to 

dSo = _8 tt% 

dH 3 ff 4 ’ 1 J 

where again, prime stands for the derivative with respect to H. We should note that the results of either considering 
cosmological constant model of the DE or the initial and ultimate thermodynamic equilibrium states of Lima et al. 
model 0 ( or briefly, dSn = 0) are obtainable by inserting wd = -1 and Q = 0 simultaneously. Since Eq. (1551) 
implies 


we use Eq. m to get 


dS D dS° D dSjj dS 2 , dH 0 
~dH ~ + dHo + dHo Zff’ 


dS 2 D __ 16n 2 QH _ dS% _ dS^_ 

dH 0 ~ 9(l+u] m )H 6 (l-n D )H 0 dH 0 dH 0 ’ 


(43) 


(44) 


where, and are introduced in Eqs. (1551) and (1551) . respectively. We have also used (fl5ll and = ji- 

to derive this equation, while H and Hq are evaluated by using Eq. m for the interacting and non-interacting 
cases, respectively. Therefore, we find an expression for the thermal fluctuations which are due to the interaction 
Q = Pd- Since it seems that the term has insignificant effects with respect to the term in the gravitational 
and cosmological setups [78| , we disregard it and inserting Eq. © into to get 


Now, using ©, we reach 


P' D ,rr 3 {dS^+dSb) 

W dH = - 8 ^- 


(45) 


4^+ATIJ H2{n ~ 2) = -^(S° d + s1 d +C), (46) 

2 n 

where C is an integration constant. Additionally, due to this fact that entropy is not an absolute quantity, C can 
be set to zero without lose of generality. S^ is evaluated in Eq. (1551) . while S° D can be obtained by integrating from 
Eq. (l25l) . This is a relation for the Hubble parameter, up to the first order fluctuations, due to the interaction Q. 
Loosely speaking, based on the first order terms of the thermal fluctuations which are due to the interaction between 
the DE candidate and the other parts of cosmos, we find the mutual relation between H and Hq in the model proposed 
in Ref. [22|. It should be noted that Eqs. (14411461) are only valid when uj d = — 1 , Q = p D while po is explained by 
Eq. ©. One can also find a general relation between an un-known Q and thermal fluctuations in the models with 
arbitrary energy density for the DE, by using Eqs. (1551) . (1501) . (1551) and (1551) . 


IV. THERMODYNAMICAL DESCRIPTION OF GRAVITATIONALLY PARTICLE PRODUCTION 

INDUCED PROCESS 


Here, we focus on the LBC model proposed in Ref. [H,|2^]. In this model, like the previous model the universe 
expansion is began from an unstable initial de-Sitter spacetime and follows the radiation and matter dominated era 
in continue. Finally, the universe expansion will reach to a perpetual de-Sitter phase which is in agreement with the 
thermodynamical equilibrium condition [dll . [25^ . In this model, there is a particle production due to the gravitational 
effects leading to an additional pressure to the Friedman equations as: 



(47) 


and 



a 


8t r(p + pc), 


(48) 
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where we considered a flat background [H, [25|. p is the energy density of the dominated prefect fluid, such as the 
radiation, confined to the apparent horizon of the FRW universe whiles, p is the corresponding pressure. The additional 
pressure (pc) plays the key role in the change of the expansion phase, and depends on the particle production rates. 
In a simple approach, considering an adiabatic process for the p artic le creation, the entropy per particle is constant 
whenever the total entropy is not [95] . Therefore, one reachs [2M03 


PC = - 


p{ 1 + w)T 
3H 


(49) 


whiles T is the particle production rate with dimension of (time) 1 
implies 


23|. 


the Energy-momentum conservation law 


p + 3Hp(l + u) = -Q, (50) 

where Q = 3Hpc [H, [25|. In the absence of this particle production (Q = 0), the entropy of the fluid with density p, 
pressure p and the state parameter ui reads as 

dS° = —-^3 (1. + 3u}°) dH o, (51) 

where we have followed the recipe of section II. Indeed, one can use Eqs. (l24l) and (fl8l) . substitute p° = 0, and finally 
replace pn with p together with using (ED to get this equation. We should note that since for the cosmological 
constant model p = 0 leading to dH 0 = 0, this equation covers the result of considering the cosmological constant 
model. We used the subscript/superscript (0) to indicate that there is no particle production in this situation. In 
addition, due to the universe expansion, the densities of confined fluids, including radiation and etc., are diluted. This 
provides a suitable situation in which the effects of this particle production process overcomes those of the other fluids 
leading to change the expansion phase. Bearing Eq. (1471) in mind, in order to evaluate the effects of this pressure on 
the entropy of the dominated fluid, we use Eq. (l38l) while p = 0 to get 

dS _ 7r(l + w) / 187i\ff 4 (l + w) — 16ttQH ^ /co ^ 

dH ~ H 3 ( 3i7 4 (l + w) 2 ’ 

where H is the Hubble parameter when the effects of the particle production is considered. In order to calculate the 
corresponding thermal fluctuations due to this pressure, by using Eq. we get 


CT 2 = 


3(1 + 3w°) 
4 


(53) 


leading to dS 1 = 0 for ui° = cons. Since this result is independent of pc , we can conclude that the gravitationally 
induced particle production processes induce weak fluctuations into the thermodynamic properties of the cosmos 
supported by a prefect fluid with constant state parameter. It explains that_why these processes do not disturb the 
current thermodynamic equilibrium state of the cosmos investigated in [24l . |25| . In order to find an expression for S 2 , 
one can insert Eqs. (ED and (15^1) into (EH1) as 


dS 2 dS dH dS 0 
dH 0 ~ dH dH 0 ~ dH a ’ 

which yields 


(54) 


dS 2 _ it(1+uj) ,18nH 2 + 6T ^ dH , tt m , 0 
dH 0 H 3 3U 2 (1 + w) 2 } dH 0 + H 3[+ 

for models obeying (l49l) . It is easy to show that 


Ho = -^o 2 (l + w°). 


(55) 


(56) 


where pc = 0, and 


H = -\h 2 ( 1 + w)[1 



(57) 
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while pc 7 ^ 0 and obeys (l49l) . In deriving these equations (Raychaudhuri equations) we have used Eqs. (l47l) and (l50l) . 
Now, since 


dH _ H 
dH 0 ~ H 0 ’ 


(58) 


we get 

dS 2 7r (1 + w ) 2 187rif 2 -f 6T r (1 + 3w°) 

dHo = + u °) ( 3H 2 (1 + w ) 2 ~~ ^ ~3LT + W a 


(59) 


Therefore, we find an expression for the fluctuations of the dominated fluid entropy which are due to the gravitationally 
induced particle production process. Finally, we should note that since the thermal fluctuation theory allows to 
evaluate S 2 up to the desired order [ill, one can use this equation to find a general relation for the particle production 
rate T, meeting Eq. HMD , in a compatible way with thermal fluctuation theory. Indeed, our study is available whenever 
the state parameter of the dominated prefect fluid be constant. Loosely speaking, in cases with non-constant state 
parameter dS 1 ^ 0 and one should use fi3l) . just the same as the section III, to find an expression for the thermal 
fluctuations of system up to the desired order. Such situations are appearing in transition eras whiles, the dominated 
prefect fluid is slowly replaced by the new one and therefore, the state parameter uj is not a constant quantity. We 
should note again that since one can consider a dominated prefect fluid with constant state parameter ui during an 
expansion era, between the two transition eras [l|, dS 1 = 0 is valid in this situation meaning that the gravitationally 
induced particle production process induces weak fluctuations into the system. These fluctuations are of the second 
and higher orders which may explain that why this model preserves its final thermodynamic equilibrium state (2~fl . [25j | . 


V. GRAVITATIONALLY PARTICLE PRODUCTION MAY MODIFY THE HORIZON ENTROPY 


Here, we are going to study the entropy of apparent horizon of the flat FRW universe, filled by a prefect fluid with 
state parameter u> = %, which is also affected by an additional pressure pc due to the gravitationally induced particle 
production process J23|. In fact, since it is shown that the DE candidate may add an additional term to the Bekenstein 
entropy of the horizon ITol - IZT^ , we are going to investigate the probable effects of the gravitationally induced particle 
production process, as the DE candidate, on the entropy of the apparent horizon. Moreover, since we work on the 
apparent horizon, we use the Hayward-Kodama temperature which is in fact an obvious generalization of the Black 
Holes temperature to the apparent horizon of the FRW universe |2l], Pipt - PHil . From the Friedmann equation m 
and m we get 

o _ 

2 HdH = —dp, (60) 

o 

and 


dp = —3 H (p + p)dt — Qdt, 


(61) 


respectively. Combining these equations together to reach 

O _ 

2HdH -\—— Qdt = —8irH(p + p)dt. (62) 

3 

Bearing the Hayward-Kodama temperature (T = — ^(1 + t^t)) in mind [40l - l43j . we can rewrite this equation as 

O _ 

T{-2HdH - —Qdt) = -4 H 2 {p + p)dt - 2{p + p)dH. (63) 

o 

Since E = pV and dV = —jj^dH, we get dE = pH~ i dH — AttH~ 2 (p + p)dt, whiles E is the associated energy 
due to the source p confined to the apparent horizon. By inserting these relations into Eq. (16311 . we obtain 

T(—2HdH - Qdt) = - dE + 2(p — p)dH, (64) 

3 7T 


leading to 


T( 


-f> 


8tt 2 
3 H 4 


Qdt) = dE — WdV, 


(65) 
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where W = is the work density plli-H fl(l |47| . Comparing this equation with the first law of thermodynamics 
in cosmological setups ( TdS A = dE — WdV), one gets 


which yields 


dS A = ~ 


s a=j+* 


( 66 ) 


(67) 


where B is the integration constant, Q = —p{ 1 + w)r, and we have used the Friedmann equation (1471) to obtain the 
second term of RHS of this equation. Using (RJTl) together with (1471) , one can rewrite the second term of RHS of this 
equation as 


(1 


H 2 


w)r, 3 

; dt = 


Tdp 


8 J p 2 (T-y/mrp)' 


( 68 ) 


Therefore, based on this equation, the DE candidate (gravitationally induced particle production process), may add 
an additional term to the horizon entropy which leads to modify the Bekenstein limit ( Sa = ^). Previously, it has 
been shown that the dynamic candidates of DE, such as the ghost dark energy model and its generalization, may 
modify the horizon entropy [l9l - l2ll |. Moreover, it is shown that any mutual interaction between the cosmos sectors 


may lead to modify the Bekenstein limit of the horizon entropy 21 . Therefore, the additional term appearing in 
this equation is in line with the mentioned attempts (Tol-lul. It is also useful to note that in the absence of a mutual 
interaction between the gravitational and baryonic fields (r = 0), the results of the previous works are recovered as 
a desired expectance [2l|, I40(-|47| . 


VI. SUMMARY AND CONCLUDING REMARKS 


In this study, we have investigated the thermodynamics of the cosmoses with the initial and final de-Sitter space- 
times. With this aim, we have showed that the state parameter of varying vacuum, as the DE candidate with varying 
energy density, cannot be equal to —1, whenever there is no mutual interaction between the fluids supporting the 
FRW background except whiles the universe is in the de-Sitter stages with H = 0 meaning that p = 0. Thereinafter, 
we have derived an expression for the entropy changes of the introduced DE model. In addition, we have payed our 
attention to a universe in which the cosmos sectors interact with each other, and we got an entropy change relation 
due to this interaction. In continue, by focusing on the model with wd = —1, we have found its corresponding 
thermodynamical interpretation. In fact, our study shows that the mutual interaction between the varyin g va cuum 
and the other parts of cosmos which disturbs the initial de-Sitter phase of the universe (whiles Ho = 0) |22j], may 
also disturb the thermodynamic equilibrium of initial de-Sitter spacetime (where dSp = d 2 S = 0), and leads to 
dSo 7^ 0. Thereinafter, since this interaction brings the ultimate de-Sitter spacetime (where again H = H 0 = 0) (22] 
which satisfies again the thermodynamic equilibrium conditions [24j], pu = 0 and thus Q = 0. The latter means that 
dSn = dSp and thus, the thermodynamic equilibrium conditions are marginally satisfied in the ultimate de-Sitter 
spacetime ( dS^ = d 2 Sp = 0) which is in agreement with J24|. Therefore, it seems that this mutual interaction 
between the initial and final de-Sitter eras may leave some thermal fluctuations into the system. Finally, by bearing 
the thermal fluctuation theory in mind, we got a relation between the thermal fluctuations of the DE candidate and 
the mutual interaction between the cosmos sectors. Moreover, we have pointed to the models in which there is a 
particle production process induced by the gravity. We showed that such process may be considered as an interacting 
model whiles this interaction comes from the gravitationally induced particle process. In fact, our resolution shows 
that such particle production, coming from non-equilibrium thermodynamic analysis, inspires a weak fluctuation to 
the thermodynamical properties of the model, including the entropy and temperature, and thus the spacetime fea¬ 
tures such as its apparent horizon. We have found out a relation between the rate of the particle production and 
these fluctuations which are second order onwards. The latter may explain why in this model the thermodynamic 
equilibrium state of cosmos does not change due to the gravitationally induced particle production process [2J, |25[. 
In addition, we have tried to establish the first law of thermodynamics on the apparent horizon and therefore, we 
found that the gravitationally induced particle production process, as the DE candidate, may add an additional term 
to the Bekenstein entropy of the apparent horizon. The latter is in line with the previous works claiming that the DE 
candidate and its interaction with the other parts of cosmos may change the horizon entropy GT" 
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